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ABSTRACT 

Mid  ultraviolet  Cameron  band  emission  from  carbon  monoxide  is  seen  in  plumes  of  Space  Shuttle  Orbiter  engine 
burns  in  low  earth  orbit.  The  observed  emission  has  been  attributed  to  chemiluminescence  from  two  and  three  step 
chemistry  of  a  minor  amount  of  methane  in  the  plume  with  atmospheric  atomic  oxygen.  DSMC  modeling  has  played 
an  important  role  in  determining  the  mechanism,  but  standard  DSMC  methods  show  significant  discrepancies  in  the 
size  and  shape  of  the  radiance.  The  differences  have  been  traced  to  a  need  to  extend  the  validity  of  scattering 
treatment  to  hyperthermal  (Erei  >  1  eV)  collision  energies.  Variable  Hard  Sphere  (VHS)  and  Variable  Soft  Sphere 
(VSS)  scattering  treatments  have  been  based  on  inverse-power-law  repulsive  inter-particle  potential  energy 
functions.  The  legacy  of  such  functionality  comes  from  the  ability  of  the  Lenard-Jones  formalism  to  produce  a 
realistic  potential  well  and  the  convenience  in  mathematical  treatments.  It  is  known,  however,  that  such  potentials 
predict  interactions  that  are  too  stiff  at  hyperthermal  energies  when  fit  to  thermal  data.  Interaction  potentials  which 
have  been  characterized  at  hyperthermal  energies  indicate  that  a  Born-Mayer  exponential  formalism  is  generally 
valid  to  interaction  energies  up  to  about  50  eV.  Born-Mayer  potential  formalism  has  been  introduced  into  DSMC 
treatment  to  extend  its  validity  into  the  hyperthermal  regime  to  define  Extended  Variable  Hard  Sphere  (EVHS)  and 
Extended  Variable  Soft  Sphere  (EVSS)  treatments.  The  addition  can  be  introduced  without  computational  penalty. 
Comparisons  between  VHS  and  EVHS  modeling  at  thermal  energies  and  EVHS  with  hyperthermal  plume  radiance 
data  demonstrate  the  validity  of  the  modification. 

1.  INTRODUCTION  AND  BACKGROUND 

Atomic  and  molecular  interactions  at  thermal  energies  are  effectively  simulated  by  the  Lenard-Jones  6-12  potential: 

V(r)  =  4£[(rjr)n-(rjrf\,  (1) 

where  V(r)  is  potential  energy,  6  is  the  well  depth,  ra  is  the  distance  where  the  potential  crosses  the  zero  energy  axis 
and  r  is  the  inter-nuclear  separation.  This  inverse-power-law  potential  includes  long  range  forces,  which  accurately 
simulate  the  potential  from  the  attractive  r 6  London  dispersion  forces.  The  shorter  range  electronic  repulsion  is 
treated  approximately  by  the  r  12  term  to  yield  a  potential  well  that  is  in  reasonable  agreement  in  depth  shape  and 
distance  of  the  well  bottom  with  those  determined  by  inversion  of  data  from  molecular  beam  scattering  experiments. 
While  the  inverse-power  functionality  for  the  repulsion  does  not  have  a  theoretical  basis  aside  from  producing 
monatonic  repulsion  it  produces  a  repulsive  wall  that  that  is  generally  reasonably  accurate  up  to  collision  energies  of 
about  0.1  to  1  eV,  covering  energies  relevant  to  the  upper  limit  to  thermal  temperatures  (-4000K). 

Such  analytical  potentials  extend  to  unlimited  distances  making  their  direct  use  in  treating  the  dynamic  interaction  of 
large  numbers  of  interacting  particles  computationally  intractable.  This  well-known  difficulty  is  managed  by  treating 
the  particles  as  Variable  Hard  Spheres  (VHS),  discrete  hard  spheres  that  vary  in  size  depending  on  the  collision 
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energy,1  the  scattering  being  isotropic  in  a  reference  frame  centered  at  the  collision  center  of  mass.  The  parameters 
used  in  describing  hard  sphere  collisions  are  shown  in  Figure  1 . 


Figure  1.  Hard  sphere  collision  between  particles  i  and  j,  shown  in  the  center  of  mass  frame.  The  collision 
cross  section,  a, ;  simply  equals  nd2,  where  d  is  the  sum  of  the  radii  of  the  colliding  particles.  X  is  the  center  of 
mass  frame  scattering  angle  which  is  equal  to  2  cos ml(b/d),  where  b  is  the  center  of  mass  impact  parameter. 

VHS  model  discrete  particle  sizes  are  defined  through  theoretical  derivation  of  viscosity  from  the  analytical 
potential.  An  important  parameter  for  determining  VHS  particle  sizes  in  a  collision  is  the  viscosity  coefficient,  co, 
which  is  linked  to  the  repulsive  exponent,  n,  of  the  inverse  power  law  potential  through  the  relation  a>  =  (n  +  2)12 n. 
The  other  parameters  which  characterize  VHS  particle  sizes  are  a  reference  diameter,  dK f,  associated  with  a 
reference  temperature,  Tref-  Reference  temperatures  are  also  identified  with  reference  relative  collision  energies 
through  the  relation  Eref  =  3kTre/2  and  reference  relative  collision  velocities  through  their  kinetic  energy 
relationship.  DSMC  including  VHS  is  effective  at  reproducing  temperature  dependence  of  viscosity  and  other 
macroscopic  observables,  such  as  plume  radiance.  Such  measurable  observables  have  been  found  to  be  remarkably 
insensitive  to  the  somewhat  unrealistic  isotropic  scattering  distribution. 

Discrepancies  are  found,  however,  in  the  ability  of  pure  isotropic  scattering  to  simultaneously  reproduce  viscosity 
and  diffusion  behavior  in  DSMC  treatment.  The  Variable  Soft  Sphere2  (VSS)  model  introduces  a  bias  to  the 
isotropic  scattering  but  maintains  a  discrete  cross  section,  allowing  viscosity  and  diffusion  behavior  to  be  accurately 
treated  simultaneously  by  DSMC. 

In  applying  DSMC  techniques  to  simulate  high  altitude  plume  radiance,  resulting  from  interaction  of  the  plume  with 
the  residual  atmosphere,  the  size  of  the  observed  radiant  field  is  larger  than  that  predicted  for  spacecraft  moving  at 
orbital  velocities.3  These  differences  were  first  noticed  in  analyzing  Cameron  band  radiance  in  plumes  from  the 
Midcourse  Space  Experiment  (MSX)  Shuttle  Plume  Observation  experiment,  described  in  reference  3.  The  Cameron 
band  system  of  CO(a  3nr  — >  X  *Z+)  is  in  the  mid  ultraviolet,  principally  between  190  to  250  nm,  and  comes  from  a 
metastable  state  that  is  not  easily  seen  in  the  laboratory.  This  emission  in  shuttle  orbiter  plumes  is  found  to  be 
dominated  by  two  and  three  step  chemistry  of  a  minor  amount  of  methane  in  the  engine  exhaust  reacting  with 
atmospheric  atomic  oxygen.  Figure  2  shows  the  MSX  satellite  image  of  steady-state  Cameron  band  radiance  for  the 
“High  ram  burn'’  reported  in  reference  3  compared  to  DSMC  model  calculations  performed  with  the  Transitional 
and  Rarefied  Axisymmetric  Monte  Carlo  Plume  (TRAMP)  code.  TRAMP  is  an  axisymmetric  DSMC  code 
developed  in  parallel  with  the  full-three  dimensional  SOCRATES  code,4  with  the  advantage  that  it  converges 
approximately  10  times  faster  for  axisymmetric  cases.  The  dimension  of  the  observed  plume  Cameron  band 
radiance  is  roughly  30%  larger  than  the  DSMC  model  predictions  for  retrofire  (into  the  ram)  burns.  Such  a 
discrepancy  is  an  indication  that  the  mean-free-path  of  interacting  species  is  not  being  treated  correctly. 


Upwind  Biased  Distance(km) 

Figure  2.  a)  Cameron  band  radiance  from  a  ram  burn  of  shuttle  orbiter  engine  at  298  km  altitude,  b)  DSMC 
simulation  of  Cameron  band  radiance,  showing  smaller  dimension  for  contour  of  a  given  fraction  of  the  peak 
radiance.  Contours  are  labeled  in  fractions  of  peak  radiance.  Atmospheric  wind  is  entering  from  left  at  46°  to 
the  view  direction.  The  cross  shows  the  location  of  the  shuttle. 


The  mean  free  path  is  a  function  of  the  particle  densities  and  the  collision  cross  sections.  Atmospheric  densities  are 
determined  using  the  MSIS  model,5  which  is  described  as  having  an  uncertainty  of  +15%.  Consistent  discrepancies 
greater  than  15%  for  multiple  plume  observations  were  seen,  so  it  is  understood  that  the  discrepancies  are  dominated 
by  errors  in  the  collision  cross  sections.  The  source  of  this  discrepancy  becomes  apparent  when  the  treatment  of 
scattering  cross  sections  by  DSMC  is  examined. 


The  usual  DSMC  scattering  treatment,  which  has  been  based  on  atomic  and  molecular  interactions  at  thermal 
energies,  deviates  significantly  from  scattering  that  occurs  at  the  hyperthermal  collision  energies  (>1  eV)  that  are 
experienced  at  orbital  velocities.  The  predicted  hyperthermal  energy  collision  repulsive  interaction  from  the  inverse- 
power  law  interaction  potentials  on  which  DSMC  scattering  has  been  based,  varies  significantly  from  the  more 
exponential-like  behavior  that  governs  the  actual  interaction.  The  importance  of  including  an  exponential 
functionality  to  the  size  of  particles  at  hyperthermal  energies  has  been  recognized  earlier  in  DSMC  treatment  of 
sputtering.6  7  Figure  2  shows  the  difference  between  the  commonly  used  Lennard-Jones  6-12  and  Born-Mayer 
exponential  repulsive  potentials  fit  at  thermal  and  hyperthermal  energies,  respectively.  While  the  Lennard-Jones 
functionality  can  fit  the  experimentally  determined  potential  at  thermal  energies,  but  the  r’12  repulsive  contribution 
deviates  strongly  at  hyperthermal  energies.  The  Born-Mayer8  exponential  repulsive  form. 


V(r)  =  A  exp (-r/B  ); 


(2) 


where  A  and  B  are  empirically  determined  range  and  softness  parameters,  respectively,  doesn’t  produce  a  well,  but 
can  fit  the  measured  potential  over  a  broad  range  of  hyperthermal  energies.  At  hyperthermal  collision  energies  the 
potential  well  will  have  little  effect  on  a  collision  outcome.  It  is  seen  in  Figure  3  that  the  Lennard-Jones  potential  fit 
to  the  He-He  potential  well  produces  much  larger  collision  turning  points  at  hyperthermal  energies  than  the 
measured  potential. 


He-He  Potential  at  Hyperthermal  Energies 


a)  b) 

He-He  Potential  at  Thermal  Energies 


Figure  3.  Lennard-Jones  6-12  and  Born-Mayer  exponential  repulsive  potential  forms  compared  to  measured 
He-He  potential  at  a)  thermal  and  b)  hyperthermal  energies. 

In  general,  an  inverse-power-law  potential  derived  to  match  the  energy  dependence  at  some  energy  will  consistently 
over  predict  representative  hard  sphere  diameters  by  an  increasing  error  as  the  collision  energy  increases. 

The  method  chosen  for  defining  variable  hard  sphere  diameters  derived  from  Born-Mayer  potential  functionality 
involves  matching  the  velocity  dependence  at  a  reference  energy  (or  a  reference  velocity  or  reference  temperature, 
as  mentioned  earlier).  The  new  diameters  and  cross  sections  are  referred  to  here  as  Extended  Variable  Hard  Sphere 
(EVHS)  diameters  and  cross  sections.  EVHS  diameters  can  also  be  used  in  combination  with  VSS  scattering 
distributions,  effectively  defining  an  Extended  Variable  Soft  Sphere  (EVSS)  applicable  at  hyperthermal  energies. 
For  valid  EVHS  and  EVSS  cross  sections,  however,  it  is  important  to  use  VHS  cross  section  parameters  from  a 
reference  energy  that  is  high  enough  that  the  interaction  is  dominated  by  exponential  repulsion  and  not  influenced 
appreciably  by  long  range  attraction.  A  prescription  for  deriving  EVHS  model  parameters  from  VHS  functionality 
and  incorporating  EVHS  cross  sections  into  DSMC  has  been  published  earlier.9  Typical  cross  section  ratios  derived 
from  EVHS  and  VHS,  matched  at  a  relative  collision  velocity  of  2.5  km/s  are  shown  in  Figure  4. 


Relative  Collision  Velocity  (km/s) 

Figure  4.  Velocity  dependence  of  typical  EVHS/VHS  cross  section  ratio  where  the  cross  sections  have  been 
matched  at  a  relative  collision  velocity  of  2.5  km/s. 


This  paper  includes  derivation  of  viscosity  and  collision  frequency  formulas  based  on  EVHS  Born-Mayer  potential 
functionality  and  verification  through  comparison  of  viscosity  and  shock  density  profiles  for  VHS  and  EVHS 
models.  Validation  at  hyperthermal  energies  is  performed  through  applying  EVHS  to  modeling  the  Cameron  band 
emission  data  shown  in  Figure  2a. 


2.  VISCOSITY  AND  COLLISION  FREQUENCY  OF  EVHS  MODEL 


For  the  numerical  results  obtained  with  the  EVHS  model  to  agree  with  experimental  data  one  needs  to  match 
viscosity  of  simulated  gas  with  that  of  the  real  one.  To  this  end,  let  us  derive  the  temperature  dependent  expression 
for  the  viscosity  of  an  EVHS  gas.  The  1st  order  approximation  to  the  coefficient  of  viscosity  in  a  single  species  gas 
may  be  written10  as 
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where  m  is  molecular  mass,  T  is  gas  temperature,  g  is  the  relative  collision  velocity,  and  ay  is  the  viscosity  cross 
section.  For  EVHS,  similar  to  VHS,  one  can  write 
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where  d  is  the  collision  diameter  that  for  EVHS  generally  depends  on  the  relative  translational  energy  Etn  since 
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Here,  A  and  B  are  the  parameters  of  the  EVHS  potential.  Substituting  Eq  (4)  into  Eq  (3),  and  conducting  change  of 
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where  C  = - .  The  integration  of  this  expression  results  in  an  analytic  form  for  the  viscosity  coefficient, 
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where  y  =  0.5772157.  Similarly,  the  expression  for  the  collision  frequency  of  EVHS  molecules  may  be  obtained  as 
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3.  NORMAL  SHOCK  WAVE 

The  first  application  of  the  EVHS  potential  presented  in  this  paper  is  the  structure  of  a  normal  shock  wave.  The 
availability  of  experimental  data11  on  shock  wave  profiles  in  argon,  as  well  as  the  measured  viscosity  data12  for 
argon  in  a  wide  range  of  temperatures,  allows  model  validation  to  be  performed.  In  addition  to  the  EVHS  model,  the 
standard  VHS  model  has  also  been  used.  Prior  to  the  DSMC  runs,  the  viscosity  data  were  analyzed,  and  EVHS  and 
VHS  parameters  were  chosen  so  that  the  DSMC  viscosity  agrees  reasonably  with  the  measured  values.  The  EVHS 
parameters  used  are  A  =  1.023xl0"l4J,  B  =  0.253  A  at  a  reference  temperature  of  273K.  For  the  VHS  model,  dref  = 
3.73  A  and  a  =  0.67.  As  was  shown  in  reference  9,  the  parameters  of  the  EVHS  model  A  and  B  may  be  related  to  the 


VHS  parameters  dref  and  co,  if  one  assumes  (£’/r)EVHs=(£/r)EVHs  and  (£’/,)'EVHs=(£fr)'EVHs  at  some  reference 

2  ^  ^ 

temperature.  In  this  case,  A=e  kTrej  (r  (2.5  —  ft)))  °5 ,  B  = - - - dref. 

The  above  EVHS  and  VHS  parameters  were  selected  to  provide  a  good  approximation  for  the  viscosity  of  the  EVHS 
and  VHS  models  as  compared  to  the  measured  profiles,  shown  in  Fig.  5.  Note  that  there  is  some  difference  at  lower 
temperatures  (below  about  500  K),  with  the  models  producing  somewhat  higher  viscosity  than  that  observed  in  the 
experiment,  but  the  difference  becomes  negligible  at  higher  temperatures.  It  is  recommended  that,  where  possible, 
EVHS  parameters  be  selected  by  direct  fitting  to  viscosity  in  the  highest  available  temperature  range.  In  some  cases, 
however,  it  may  be  that  some  assumed  VHS  parameters  are  available  even  though  the  fundamental  transport 
properties  are  not.  In  fact,  this  is  the  case  for  the  SOCRATES  example  shown  below.  In  this  case,  the  relations 
above  provide  the  means  to  obtain  EVHS  parameters  from  given  VHS  parameters.  We  note  strongly  that  the  two 
sets  of  cross  sections  will  cross  at  the  specified  value  of  Eref  (sometimes  specified  as  Tn.j)  and  diverge  at  higher 
energies,  so  the  conversion  should  be  done  at  the  highest  possible  value  of  Eref.  A  low  reference  value,  such  as  Tref~ 
273K,  will  match  the  cross  sections  in  a  thermal  regime  which  may  be  substantially  influenced  by  attractive  forces 
and  will  reduce  the  validity  of  the  EVHS  assumptions.  For  example,  the  VHS  parameters  recommended  for  Ar  in 
reference  10:  Tref=  273K,  c/„,;  =  4.17,  and  co  =  0.81,  begin  to  deviate  from  the  measured  viscosity  in  the  range 
T>600K.  If  the  VHS  parameters  recommended  above  for  the  M  =  3.8  shock  are  converted  to  their  high-temperature 
equivalent  dref=  3.09  and  co=  0.67  at  Tref  =  2500  K,  then  the  successive  EVHS  conversion  of  the  latter  would  predict 
A  =  6.67e-15  .1,  B  =  0.2634.  These  produce  similarly  good  agreement  with  known  viscosity  in  the  600-3000K  range 
as  the  A  and  B  values  obtained  from  direct  fitting  of  the  viscosity. 


Figure  5.  Comparison  of  DSMC  and  measured  viscosity  in  argon. 

Comparison  of  the  computed  and  measured  normalized  density  profiles  for  M  =  3.8  is  given  in  Fig.  6.  In  this  plot, 
similar  to  those  in  reference  11,  the  distance  variable  is  normalized  by  the  mean  free  path  in  front  of  the  shock  wave, 

A  =  1.098mm  for  T  =  300  K  and  p  =  50  mTorr;  the  gas  density  is  normalized  as  — — ^  .  It  is  clearly  seen  that 

Poo  -  P-0 O 

the  computed  density  values  are  very  close  to  the  measured  one.  There  is  a  small  difference  in  the  shallow  part  of 
the  shock  wave,  with  the  computed  profiles  appearing  to  be  more  viscous,  which  is  attributed  to  higher  DSMC 
viscosity  at  lower  temperatures  (see  Fig.  5). 


XIX 


Figure  6.  Computed  and  measured  density  profiles  inside  a  shock  wave. 


4.  VALIDATION  AT  HYPERTHERMAL  ENERGIES 

The  data  used  for  validation  ofEVHS  at  hyperthermal  energies  is  the  “High  ram  bun”  shown  in  Figure  2a.  The  High 
ram  burn  was  a  37  second  retrofire  burn  of  a  single  6000  lbf  thrust  space  shuttle  Orbital  Maneuvering  System  engine 
at  an  altitude  of  298  km  in  the  earth’s  umbra.  Atmospheric  conditions  determined  by  MSIS  were:  a  total  number 
density  of  2.78  x  108  cm’3,  a  temperature  of  717  K  and  a  mole  fraction  composition  of  O/N2/O2::0. 904/0. 093/0.003. 
More  characterizing  parameters  are  given  in  greater  detail  in  reference  3. 

In  the  original  analysis  Cameron  band  radiance  was  determined  through  agreement  with  DSMC  modeling  to  result 
from  two  and  three  step  chemistry  from  a  minor  amount  (1%  mole  fraction)  of  methane  in  the  plume: 


CH4  +  O  ->■  CH3  +  OH  (Rl) 

CH3  +  O  ->  CO(a)  +  H2  +  H  (R2) 

CH3  +  O  ->■  CH2  +  OH  (R3) 

CH2  +  O  ->■  CO(fl)  +  H2  (R4) 

CO(a)  —>  CO  +  hv  (R5) 


Reactions  (R2)  and  (R4)  are  previously  unknown  chemistry  which  were  assigned  activation  energies  of  100  and  50 
kcal/mole,  respectively,  to  obtain  agreement  with  plume  shapes.  Spectator  stripping  product  scattering  distributions 
were  uses  for  Reactions  (Rl)  and  (R3).  That  distribution  for  Reaction  (Rl)  was  found  to  be  important  in  reproducing 
the  observed  plume  shape,  and  subsequently  has  been  verified  in  hyperthermal  molecular  beam  laboratory 
experiments.13  The  remaining  discrepancies  in  plume  size  and  shape  predicted  by  plume  modeling,  shown  in  Figure 
2,  were  the  inspiration  for  postulating  the  importance  of  the  EVHS  treatment.  The  axisymmetric  DSMC  TRAMP 
code  was  modified  to  use  EVHS  cross  sections.  VHS  reference  cross  sections  at  reference  velocities  ranging  from  2 
to  8  km/s  were  used  to  generate  EVHS  characterizing  parameters.  Figure  7  shows  the  comparison  of  the  data  with 
the  EVHS  computed  plume,  with  all  the  remaining  input  parameters  remaining  the  same,  showing  a  significant 
improvement  in  agreement  in  the  size  of  the  contours. 


Figure  7.  Comparison  of  observed  Cameron  Band  radiance  with  EVHS  modeled  radiance. 

The  TRAMP  code  uses  a  global  value  of  the  viscosity  coefficient  co.  The  value  that  had  been  used  predominantly  in 
past  studies,  including  the  cases  shown  in  Figures  2b  and  7,  was  co  =  0.75.  Before  including  the  EVHS  model  into 
TRAMP  an  attempt  to  increase  the  mean  free  path  for  molecules  traveling  at  hyperthermal  velocities  was  done  by 
increasing  co  to  0.9  without  changing  the  species  reference  diameters  and  energies.  That  change  indeed  increased  the 
size  of  the  radiance  contours  but  also,  interestingly,  modified  the  shape  of  the  0.5  (of  peak  radiance)  contour  to  more 
closely  match  its  “egg”  shape  in  the  data.  It  was  anticipated  that  the  EVHS  addition  might  accomplish  the  same 
improvement  in  shape.  The  comparison  in  Figure  7  shows  a  small  improvement  in  shape  but  indicates  that  there  are 
still  discrepancies  that  indicate  a  larger  co  may  relevant  to  the  collisions  governing  the  shape  of  the  radiance.  The  co 
=  0.75  value  was  determined  from  examining  viscosity  coefficients  for  stable  gases,  all  of  which  have  closed 
electronic  structures.  Radicals,  such  as  atomic  oxygen  have  unpaired  electrons  which  would  be  expected  to  produce 
softer  repulsive  potentials  in  collisions  with  closed  shell  structures.  Considering  that  the  spatial  distribution  at  the 
level  of  the  0.5  contour  is  probably  most  strongly  influenced  by  O  +  CH4  collisions  co  =  0.8  was  use  for  an  EVHS 
analysis.  The  result,  shown  in  Figure  8,  shows  an  improved  match  in  shapes  of  the  0.5  contour  and  generally 
improved  agreement  among  all  the  contours. 


Figure  6.  Comparison  of  observed  Cameron  Band  radiance  with  EVHS  modeled  radiance  with  co=  0.8.  Data 
have  been  offset  ~2  km  to  the  right,  consistent  with  the  shuttle  beacon  tracking  uncertainty. 


5.  CONCLUSIONS 


A  simple  modification  of  the  VHS  collision  model  is  described  that  will  be  useful  for  DSMC  cases  where  the  typical 
collision  energy  greatly  exceeds  the  collision  energies  typical  of  thermal  viscosity  measurements.  For  flowfields 
with  hyperthermal  collision  energies,  the  EVHS  and  EVSS  can  be  expected  to  provide  a  superior  method  of 
extrapolation  from  known  thermally-based  VHS  parameters.  Care  must  be  taken  in  applying  this  model,  since  it 
assumes  that  the  VHS  parameters  used  as  a  starting  point  reflect  a  regime  dominated  by  a  repulsive  potential.  If  the 
VHS  parameters  are  based  on  viscosity  values  that  reflect  a  significant  contribution  from  attractive  forces,  then  the 
validity  of  this  method  of  predicting  the  high-energy  collision  cross  section  is  lost. 

We  have  verified  the  ability  of  this  model  to  continue  to  reproduce  the  given  viscosity  data  near  the  reference 
temperature.  Validation  data  for  very  high  temperature  or  very  high  energy  flows  are  rare,  but  we  have  provided 
one  comparison  case  with  measurements  that  validate  the  ability  of  this  model  to  improve  the  physical  realism  of  the 
collisions  predicted  by  DSMC.  One  day,  it  may  be  that  any  viscosity-based  simplified  collision  model  will  be 
unnecessary  as  computational  chemistry  tools  may  provide  reasonable  values  for  elastic  collision  cross  sections  (as 
well  as  inelastic  cross  sections)  over  the  full  range  of  energies.  First  examples  of  this  capability  have  been 
demonstrated  by  Tokumasu  and  by  Ozawa  [refs]. 
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